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Abstract 

In this paper we explore the magnetic and orbital properties closely related to a tetragonal-orthorhombic 
structural phase transition in iron pnictides based on both two- and five-orbital Hubbard models. The 
electron-lattice coupling, which interplays with electronic interaction, is self-consistently treated. Our re- 
sults reveal that the orbital polarization stabilizes the spin density wave (SDW) order in both tetragonal and 
orthorhombic phases. However, the ferro-orbital density wave (F-ODW) only occurs in the orthorhombic 
phase rather than in the tetragonal one. Magnetic moments of Fe are small in the intermediate Coulomb in- 
teraction region for the striped antiferromangnetic phase in the realistic five orbital model. The anisotropic 
Fermi surface in the SDW/ODW orthorhombic phase is well in agreement with the recent angle-resolved 
photoemission spectroscopy experiments. These results suggest a scenario that the magnetic phase transi- 
tion is driven by the ODW order mainly arising from the electron-lattice coupUng. 

PACS numbers: 75.30.Fv,74.20.-z,71.10.-w 
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I. INTRODUCTION 



The recent discovery of superconductivity in RFeAsi_;fOr (R=La, Ce, Sm, Pr, Nd, etc.) with 
high transition temperature has attracted extensive interest yj]. The parent compound LaFeAsO 
shows strong anomalies near 150 K in resistivity, magnetic susceptibility, specific heat, etc.. It is 
suggested that the ground state is a spin-density-wave (SDW) ordered state with a striped antifer- 
romagnetic (S-AFM) configuration [|2|]. The predicted magnetic structure was then confirmed by 
a subsequent neutron diffraction experiment [|3|]. Interestingly, the neutron diffraction data indi- 
cated that a subtle structural distortion occurs first at ~ 155 K, and the SDW long-range order 
establishes at a slightly lower temperature, r5/)vi/ ~ 137 K [j]. This phenomenon also has been 
generically found in other 1111 systems of the iron pnictides, e.g., NdFeAsO {T^ and Tsdw ~ 150 
K and 141K, respectively.) |4, ^, and in the 111 systems, e.g. NaFeAs with 52 K and 41 K 1^], 
and in the 11 systems, e.g. FeTe with 87 K and 75 K [7]. However in 122 systems, the structural 
and magnetic phase transitions spontaneously happen, i.e., at 142 K for BaFe2As2 [8], and at 205 
K for SrFe2As2 ll9l-lllll. etc.. Although the structural phase transition occurs before the magnetic 
ordering is formed, it is difficult to distinguish either the instability of the electronic structure 
or the orthorhombic lattice distortion plays an essential role. Especially, unless there exists strong 
magneto-elastic coupling, it is very strange that the structural transition and the magnetic transition 
are so closely related. 

The SDW order is suppressed by electron or hole doping, and the superconductivity is triggered 
in doped iron pnictides [|2|]. The closeness of the superconducting phase to the SDW instability 
implies that the AFM fluctuations play a key role in the superconducting pairing mechanism. 
Therefore, investigating the origin of the SDW instability in the parent compound is an essential 
step to understand the microcosmic origin of the superconductivity. Up to date, two main mecha- 
nisms were proposed for addressing the SAFM ordering in iron pnictides. The first one suggested 
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1511 . which is 



that the SAFM results from the Fermi-surface nesting of itinerant electrons [|2i 
itinerant SDW character with modulation wavevector Q=(;r,0). Alternatively, it was proposed that 
the superexchange interaction mediated through the off-plane As atom dominates the spin config- 
uration formation [Il6l-l2lll. In such a local picture, the SAFM ordering arises from the competition 
between the nearest-neighbor (NN) and next nearest-neighbor (NNN) spin couplings, when the 
NNN exchange coupling becomes larger than half of the NN exchan ge i nteraction. Whether an 
itinerant picture [|22l - l24ll or a localized superexchange mec 

hanism[[ljQ,Q,Q] responsible for 
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the SDW ordering is still a hot debate for Fe-based layered systems n25\\ . 

On the other hand, it is now believed that iron pnictides are the multi-orbital systems, with Fe^"^ 
ion in a tetrahedral crystal field (CF). The energy level splittings of the five 3d orbitals are small. 



in a magnitude order of 0.1 eV [|13 



ni . Thus the splitting is sensitive to the lattice distortion. 
It is found that in the high-temperature tetragonal phase, the xz and yz orbitals are degenerate 
ll^l . While we expect such a degeneracy is removed in the low-temperature orthorhombic phase. 



Up to date, the minimal two-orbital model with half-filling 



a filling of one third (i.e., two e 



four electrons in three orbitals) [|30 



ectrons in three orbitals) 
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28h . three-orbital model with 



and a filling of two thirds (i.e., 



3311 . four-orbital model with half-filling [128 



, and full five- 



orbital model (with six electrons in five orbitals) in the two dimensional (2D) case Hi 51. 1321-13 811 and 
three dimensional (3D) one were proposed for addressing the low-T electronic. 



magnetic 



an d op tical properties of the iron pnictides. However, in most of previous works [|27 
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3611 . the orthorhombic CF splitting was neglected. Only a few authors suggested the role of the 



orthorhombic CF splitting ll4ll - l43ll . but the magnitude of the CF splitting was not self-consistently 
determined. Therefore, whether the degeneracy of the orbitals and the orthorhombic distortion 
play an essential role on the groundstate properties or not is still not very clear. 

The remove of the degeneracy of the xz/yz orbitals may be associated with the orbital ordering 
(00) in iron pnictides. In a recent angle-resolved photoemission spectroscopy (ARPES) exper- 
iment on BaFe2As2, it was found that magnetostructural transition is accompanied with orbital- 
dependent modifications in the electronic structure ll44ll . which is an evidence of the orbital polar- 
ization or 00. Especially, the infrared phonon anomaly in BaFe2As2 was also thought to be the 
consequence of the ordering of the orbital occupation li5\. In addition, the magnetic [46], resis- 
tance 1I47I-I49I1 and optical conductivity nSuA anisotropics were observed in the experiments, which 
were attributed to be the evidence of the orbital polarization. Meanwhile, the possible influence 
of the orbital polarization has also been investigated by many other experimental and theoretical 



The 00 concept based on the 



works, including local-density approximation (LDA) v33[ 
strong correlation and localized picture, which is usually applied for classical insulator materials, 
had been proposed to being responsible for the SDW order [jsil. However, the iron pnictides, be- 
ing a bad metal and in the moderate correlation regime [|53l[ . it doesn't seem to fall into this class 
material. Therefore, the orbital density wave (ODW) based on itinerant scenario can be applied on 



these moderate correlation systems [[54 



Although some theoretical works proposed the 00 scenario [|34l l5lll. to interpret the origin of 
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the orbital polarization of iron pnictides observed in experiment, other authors suggested that these 
results are mainly driven by the magnetic ordering in absence of the 00 bsll . Recently, the strong 
electron-lattice (e-1) interaction associated with the orthorhombic distortion in BaFe2As2 has been 



suggested to take responsiblility for the structural phase transition [|58|], which contradicts with 
the scenario that the 00 induced by purely electron-electron (e-e) interaction drives the structural 



phase transition [15 ill . Therefore, the origin of the orbital polarization and its consequence are still 



an open and debating question. As an altemative, we suggest recently that the ODW favors the 



existence of the SDW in iron pnictides [|57ll . In this paper, we propose that the SDW ordering 
is driven by the density wave-type 00 state in iron pnictides. We find that such an 00 phase 
is stabilized only in the presence of the orthorhombic distortion. The 00 physics in five orbital 
model is similar to the degenerate two-orbital picture due to the Jahn-Teller-type distortion in the 
itinerant background. This paper is organized as follows: a model Hamiltonian and mean-field 
approximation are described in Sec. 11; then the numerical results and discussions are presented in 
Sec. Ill; the last section is devoted to the remarks and summary. 

II. EFFECTIVE HAMILTONIAN AND ORDERING PARAMETERS 

We start with an extend multi-orbital Hubbard model Hamiltonian including both the e-e and 
e-1 interactions, 

H = Ho + Hi + He-,. (1) 
Here Hq describes the kinetic energy term, 

^0 = ^ ?yC,:^^Cj/j^-;/^n,a^ (2) 

i,j iacT 

where C^^^^ creates an electron on site i with orbital a and spin cr, t"^ is the hopping integral 
between the i site with a orbital and the j site with (5 orbital, and /i is the chemical potential 
determined by the electron filling. The electronic interaction part reads, 

i.a i i,(T 

-Jh ^{C*^^CiiiC^2iCi2^ + h.c.) + Jh ^{C^^^C'l^^Ci2iCi2^ + h.c.) (3) 
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where U(U ) denotes the intra-(inter-)orbital Coulomb repulsion interaction and Jh the Hund's 
coupling. We take U'=U-2J// throughout this paper. 
The e-1 interaction is depicted as 

He-i = ^ J]^nn^ - na^)e'^'-^' + (4) 

icr i 

where 6 is the CF level splitting of the xz and yz orbitals induced by orthorhombic distortion, and 
Qo is the ODW vector, A is a constant coefficient. Since the metallic iron pnictides do not fall into 
conventional 00 scenario for insulators, it is only a SDW-type spin ordering and ODW-type 00 
for complicated metals [Svl. In the presence of the SDW and ODW order, we adopt the following 
mean-field approximation in real space to decouple the particle-particle interaction terms in Eq. 
(3), 

< "iacr > = -[ni2 + (Tm,e'^'''^' +amae'^'''^' +cram,oe'^^'*^''^'^']. (5) 

For the a orbitals (l=zx and 2=yz) in iron pnictides, and for the electrons in /3 orbitals with 
/3 = 3 : - y^,4 '■ xy,5 : 3z^ - r^, we adopt, 

1 ^ ^ 

< > = -[nfi + (6) 

where we define the particle numbers nn = jj Y,k<T,a=i,2 < C^^a-^kacr >,n/^ = Y^ko- < C^/so-^kpa- >, 
and the SDW/ODW order parameters m, = jj Z/tcr,a=i,2 cr < C'l^^Ck-Q^aa- >, rrio = jj Zicr,a=i,2 Q' < 

^kaa-^k-Q„ao- >, m^o = jj 'Lko;a=i,2 ^'^ < ^kao-^ l^-Q^~Qoa<^ >' "^Z? = jj Tjko- < C'l^a-^k-Q,/3o- >, 

respectively. In the momentum space the decoupled eff"ective Hamiltonian can be written as 

Ho = Yj^T"''(^^Cta.Ck/S<r + Tf'"ik)Cl^^Cka.-I^Cl^Cka^^ (7) 



and 



^/ - ^ [{ACl^^CkacT + 0-BCl^^^Ck-Q,a:T (8) 



k,o- 

a=zx,yz 



+«CCl^^Ck-Q„aa + O-aDCl^^Ck-Q^-Q^aa)] 

+ J] [{EfjClp^Ckpa + o-FfjClp^Ck-Q.fia-)] + const. 



kcT 



where the coefficient A = ^nuU + ^{^nu + + n4. + ns){2U' - Jh), B = -^nisU - \{\ms + + 
Ma + ms)JH, C = \mo[U - W + Jh],D = -\mso[U - Jh], Ep = \npU + {{nn + n0{2U' - Jh), 



and Fp = -^m/jU - \{ms + mp)Jii- The constant term reads 

const. = ~'^^[2^"i2 ~ ~ ^so) + «3 + «4 + ^5 - - m4 - mj] (9) 

-//i7 [-(«i2 - W^o) + ^12(^3 + ^4 + "5) + '^3'^4 + ^3/^5 + '^4/^5] 

+ Y-^//[^(«?2 + - "^o - "^L) + «12(«3 + «4 + "5) + «3«4 + «3«5 + "4«5 

+ms(m3 + m4 + nis) + m3m4 + msms + m^ms^ 
The e-1 interaction describing the orthorhombic distortion is expressed as 

He-l = 2 y^j(C'l,xcrCk-Q,zx<T " C'ly^^Ck~Q„y,,a-) + —5^ (10) 

Minimizing the groundstate energy with respect to the energy level splitting gives rise to the 
level splitting 5 = Y^kaa-'^ < ^taa-^k-Qoaa- >■ Thus the CF splitting 5 can be obtained self- 
consistently hy 6 = -^nio/A. 



III. RESULTS AND DISCUSSIONS 



In this section, we first present our studies for the two orbital model, and then for the more 
realistic five orbital model of the quasi-two-dimensional systems, such as 1111 phase (LaOFeAs, 
etc.), within the mean-field approximation, finally extend to the three-dimensional systems in the 
122 phase (BaFe2As2, etc.). The magnetic and orbital properties closely related with a tetragonal- 
orthorhombic structural phase transition in the iron pnictides are explored. In addition, the effect 
of the inter-player coupling for BaFe2As2 is also considered in the five orbital model. 



A. Two-Orbital Model 

We firstly present our results in the case of the two-orbital model [Q]. The possible spin and 
orbital configurations with modulation wave-vectors Qv/Qo=(0,0),(0,;r),(;r,0) and (n,7r) in both the 
tetragonal structure and the orthorhombic phase are studied. Note that in the tetragonal phase, the 
orthorhombic distortion is absent, i.e. H = Hq + Hj\ for comparison, while in the orthorhombic 
phase the orthorhombic distortion is present, i.e. H = Hq + Hj + H^^i. The J//-U phase diagrams 
are obtained both for the tetragonal phase and for the orthorhombic one, as shown in Fig. 1 (a) and 
(b). In the tetragonal case, the SAFM metallic phase is a stable ground state with small magnetic 
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moment in the intermediate Coulomb interaction parameter region, and all the orbital configura- 
tions are degenerate in the SAFM phase, as seen in Fig. 1 (a). So the ground state is nearly a 
para-orbital or orbital liquid phase. However in the orthorhombic phase, the e-1 coupling breaks 
the orbital degeneracy of the xz-orbit and the yz-orbit, leading to a ferro-orbital (FO) configura- 
tion. Due to the itinerant character of the electron, the orbital polarization in the ODW state is 
small. Since the pure electronic interaction is insufficient to contribute the FO-SAFM phase, the 
Jahn-Teller-type orthorhombic distortion removing the xz- and yz-orbital degeneracy results in a 
weak FO ordering phase. While, in the large U parameter region, the 00 is destroyed under the 
strong Coulomb interaction, which favors electron occupation in the xz and yz orbitals equally. 
In addition, the FO phase absence of SDW appears in the weak Hund's coupling region, which is 
the result of the competition between the crystal field splitting and the Hund's copling. Here, the 
elastic coefficient parameter A is fixed to be 0.5/t (t is the hopping integral parameter in Ref. 924^ ) 
for the two-orbital model throughout this paper. 



0.30 



0.05 




(b) 



U/t 




FIG. 1: I//-U phase diagram of the two-orbital model in the tetragonal (a) and orthorhombic (b) phases. 
PM and SAFM denotes paramagnetic and striped antiferromagnetic iQs=(n,0)) phase, respectively. FO 
represents the ferro-orbital order with 2o-(0, 0). The solid lines with different symbols display different 
phase boundaries. The dash line and short dash line denote the the metaf-insufator border and the fow 
magnetic (LM) moment(// < 1 //g) and high magnetic (HM) moment (p> I /ib) phase border, respectively. 



Minimizing the groundstate energy, we self-consistently obtain the orthorhombic CF splitting, 
S = -jifio/A, for the orthorhombic case. In Fig. 2, the dependence of the CF splitting between 
the xz- and yz-orbit on the elastic coefficient A is plotted with the Coulomb repulsion U=4? for 
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different Hund's coupling J//. For small coefficient A, the CF splitting 6 becomes larger under the 
strong e-1 interaction. However, the large Hund's coupling obviously suppresses the CF splitting, 
since the Hund's coupling tends to distribute the electrons in two orbitals equal-weightly. This is 
also seen in Fig. 1(b), the FO phase only exists for small J//. 



0.0 
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— J„=0.0 

- -J^=0.075U 
^J„=0.1U 
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FIG. 2: Dependence of the crystal field splitting 6 in the orthorhombic case on the elastic coefficient A in 
two-orbital model, with the parameters U=4t for J//=0, 0.075U ,0.1U and 0.2U. 



It is known that the SDW, ODW and combined spin-orbital density wave order parameters, m^, 
irio and m^o, crucially depends on the Coulomb interaction. The comparison between the tetragonal 
and the orthorhombic phases are displayed in Fig. 3. Although the ordering parameters are rather 
small, an obvious orbital polarization occurs in the orthorhombic case in the intermediate Coulomb 
interaction range. As a comparison, it is absent in the tetragonal case. A rather small SAFM 
magnetic moment is also obtained in the intermediate coupling region, as seen in the shadow 
region of Fig. 3. Notice that with the increase of U, the crossover of the orbital polarizations from 
the xz-type symmetry to the yz-type one occurs at ~ 4.5t. 

To compare partial density of state (PDOS) in both tetragonal and orthorhombic cases are 
plotted in Fig. 4. It is shown that there is no orbital polarization in the non-interacting case at U=0. 
With the increase of U to a region of Ui~2.8t<U<Uc, the PDOS of the xz-orbital component is 
larger than that of the yz-orbital component near Fermi surface in orthorhombic case, indicating 
an obvious xz-orbital polarization. However, further increasing U to the region of Uc<U<U2~5.8t, 
it changes to a weak yz-orbital polarization, in agreement with Fig. 3. In the orthorhombic phase, 
there exists larger orbital polarization than that in the tetragonal phase. Obviously, although the 
electronic interaction plays a key role in the orbital polarization, the e-1 interaction finally stabilizes 




FIG. 3: Dependence of the spin, orbital and spin-orbital ordering parameters on Coulomb interaction in the 
tetragonal (6=0) and orthorhombic {6i=0) phases, with J//=0.1U. The vertical dashed dotted Une denotes the 
boundary of different orbital polarization. The shadow region represents the small magnetic moment case 
(< 1 I^b)- 

the ground state of the LaOFeAs to the FO-type SAFM phase. 

To get a further insight into the SDW and ODW states, we display the evolution of the Fermi 
surface on the Coulomb interaction in the folded Brillouin zone of SAFM phase in Fig. 5. The 
Fermi surface nesting is obviously observed in the non-interacting paramagnetic case, including 
two hole pockets at the center and two electron pockets at the comers of the Brillouin zone. Never- 
theless, with the increasing of the electronic interaction, the nesting becomes weak and the ordered 
SAFM state is stabilized. 

In addition, both the hole doping and the electron doping behaviors are investigated in this two- 
orbital model, the results similar to Ref.[27] are obtained. Comparing with the experiments, we 
find that in the hole doped case, the magnetic moment decreases with doping, and is concordant 
with the experimental results While the AFM moment increases with the increase of electron 
doping, which is completely in contradiction to the experimental observations [|2|]. On the other 
hand, the Fermi surface of the ordered state and the xz-orbital polarization character of the Fermi 



surface in the two-orbital model are not in agreement with the recent ARPES experiments [|44(|. 
To resolve these discrepancies, we have to extend the present effective two-orbital model to the 
realistic five-orbital model. 
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B. Five Orbital Model 



It is now generally believed that the more orbitals other than two orbitals are involved in the 
low-energy physics in iron pnictides. The five orbital of 3d electrons of Fe^^ in the tetrahedral 
CF should be considered to address the electronic, magnetic and orbital properties of the iron 
pnic tides. In this subsection, we explore the role of the e-1 coupling within the five-orbital model 
OSn . Notice that, unless otherwise specified, the elastic coefficient parameter A is fixed to be 1.0 
eV"' for the five-orbital model throughout this paper. Phase diagram of the five-orbital model for 
the tetragonal (a), and orthorhombic (b) phases are shown in Fig. 6. It is clearly found that the FO 
phase is in a stable ground state in the presence of the orthorhombic distortion. The intermediate 
Coulomb interaction with a weak Hund's coupling favors the SAFM phase with low magnetic 
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FIG. 5: Evolution of Fermi surface of the two-orbital models in the folded Brillouin zone (a) U^O, (b) 4t, 
(c) 5t, and (d) 6t with Jh= O.IU for the orthorhombic phase. 

moments, similar to the two-orbital model. However, the five-orbital results are obviously different 
from that of the two-orbital model in Fig. 1(a) and Fig. 1(b). With the increase of the electronic 
Coulomb interaction, the Neel AFM phase with Qs=i7T, n) is in the ground state. Moreover, in 
strong Coulomb repulsion and large Hund's coupling region, the ground state of the system is the 
FO-SAFM phase with large magnetic moment (//>1 //g) due to the Hund's coupling. This result is 
obviously different from the one without FO in the simple two orbital case, which implies the five 
orbital case involves a complicated multi-orbital correlated effect. Further increasing the Coulomb 
correlation leads the system transit to insulator. 

The dependence of the magnetization of each orbital, and the orbital occupancy, on the 
Coulomb interaction in both the tetragonal and the orthorhombic structure are obtained, the or- 
thorhombic case is shown in Fig. 7. One finds that the xz-orbital polarization is larger than that 
in the tetragonal one in a very narrow intermediate Coulomb repulsion region from U=1.2 to 1.45 
eV for 7//=0.25U. In this region, the low magnetization, 0</z<l hb, with bad metallic state is also 
obtained. The electron occupancy in the xz-orbit is always more than that in the yz-orbit as the 
Coulomb interaction increases. No orbital occupancy crossover from the xz-type to the yz-type 
symmetry is observed, which is different from the two-orbital results above. This shows that the 
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FIG. 6: J//-U phase diagram of the five-orbital model for the tetragonal (a) and orthorhombic (b) phases, 
respectively. PM, SAFM and AFM denotes paramagnetic, striped antiferromagnetic iQs=i^,0)) and Neel 
antiferromagnetic (2.v=(7r, tt)) phases, respectively. FO represents the ferro-orbital order with 2o-(0, 0). 
LM and HM denote the low magnetization (ju<l hb) and high magnetization Ou>l //b), respectively. 

ferro- xz-orbital polarization arises from the electronic interaction of the five 3d orbitals, while the 
long range ordering, i.e. 00 is contributed from the e-1 interaction. Our results are consistent with 



not only the crystal structure but also the ARPES experimental results [|44ll . where the xz-orbital 
polarization is observed near the T point. 

In the present five-orbital model, the orthorhombic CF splitting of the xz- and yz-orbit is also 
determined self-consistently. The dependence of the CF splittings on the e-1 coupling, shown in 
Fig. 8, is similar to that in two-orbital situation. On the other hand, in contrast to the two-orbital 
case, where the Hund's coupling suppresses the CF splitting, the influence of the Hund's coupling 
on the CF splitting is complicated in the five-orbital model because of the multi-orbital effects. Ei- 
ther the strong or the weak Hund's coupling J// favors the large orthorhombic CF splitting. Hence, 
it indicates that the orthorhombic distortion is associated with the multi-orbital effect. Thus, due 
to the xz-orbital polarization with the ODW of wave vector (0, 0), it suffers an orthorhombic dis- 
tortion, leading to that the lattice parameter a is larger than b, as experimentally observed. 

In Fig. 9, the PDOS are plotted for various magnetic configurations: non-magnetic, FO- 
SAFM(LM), Neel AFM without 00, and FO-SAFM(HM), respectively. In the non-interacting 
(a) and FO-SAFM (b) cases, three \.2g orbitals, xz, yz and xy components, mainly contribute to 



the FS, which agrees with the LDA calculations [|26|] and the ARPES experiments [|59|1. In non 
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U (eV) U (eV) 

FIG. 7: Dependence of total magnetic moment, orbital magnetization (a), and orbital occupancy (b) on 
Coulomb interaction in the orthorhombic structure for J//=0.25U. The shadow region denotes the small 
magnetic moment physical parameters region with jU<l jip. 
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FIG. 8: Dependence of the crystal filed splitting of xz and yz orbitals on the elastic coefficient with U=1.4 
eV for diff'erent Hund's rule coupling. Other parameters are the same with Fig.7 

magnetic case, there is no orbital polarization of xz and yz orbitals in FS. While in the small 
magnetic moment phase, the PDOS of the xz-component is larger than that of the yz-component, 
indicating the orbital polarization is mainly xz-component, as seen in Fig. 9(b). In the Neel AFM 
phase with Qo=(7r, n), the - orbital mainly contribute to the FS . On the other hand, in the FO- 



14 



SAFM(HM) phase, there is also an obvious xz-orbital polarization near the FS. Notice that due to 
the large magnetic moment and strong Coulomb interaction, the electrons are localized electrons 
rather than itinerant electrons in the SAFM states. 
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FIG. 9: Partial density of states in the paramagnetic (a), FO-S AFM(LM) (b), Neel AFM (c), and FO-SAFM 
(HM) (d) phases in the orthorhombic phase, with parameters U^O, 1.25, 1.75 and 2.25 eV, respectively. 
The Hund's coupling J//=0.25U is adopted. 



With the increase of the Coulomb interaction, the evolution of the FS is displayed in Fig. 10. 
From the PM state to the FO-SAFM state, the FS nestings are destroyed with the increasing elec- 
tronic correlations. However, in the Neel AFM states, the FS nesting remains, which leads to the 
(n, n) AFM states. With the further increasing strong Coulomb interaction, the FS nestings become 
weak. Especially, the anisotropy of the FS appears in the SDW/ODW state at the Coulomb inter- 
action U~1.25 eV, with a larger hole pockets along r-F(0,;r) in the folded SDW/ODW Brillouin 
zone corresponding to the ferro-magnetic direction and a smaller one along r-X(7r,0) correspond- 
ing to the anti-ferromagnetic direction in the SAFM state. The anisotropic character of the FS 
resembles the results of the recent ARPES experiments [|60||. The anisotropic FS is obviously 
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different from the FS observed in other ARPES experiments [^61^ and the de Haas-van Alphen ex- 
periment [62] . Furthermore, with the increase of U to 1 .4 eV, the anisotropy becomes weak, with 
the FS along Y-Y direction disappearing. Our results reveal that the anisotropic character of the 
FS and other properties are mainly due to the splitting of xz and yz orbitals in the orthorhombic 
phase. Hence the anisotropy is obviously manifested by the e-1 coupling. 
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FIG. 10: Evolution of Fermi surface in the nonmagnetic/paramagnetic (a) and (b), FO-SAFM of LM (c) 
and (d), Neel AFM (e), and FO-SAFM of LM (f) phases in the orthorhombic phase, with the Coulomb 
correlation parameters are U=0 (a), 1.0 (b), 1.25 (c), 1.4 (d), 1.75 (e), and 2.25 (f) eV. J//=0.25U is adopted. 



To compare with the two-orbital model, we also investigate the doping case in the five-orbital 
model. We find that either hole or electron doping results in the decrease of the magnetic moment, 
which is well consistent with the experimental observations [2]. It is also implied that the five- 
orbital model is a more accurate realistic model which well describes the low-energy physics in 
the iron pnictides. 
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IV. SUMMARY 



In this work, we also study the realistic five-orbital model [13911 for BaFe2As2 of more 3D sys- 
tems, with the SDW vector Qs={n,0,n), and ODW vector 2o=(0,0,0) observed experimentally. To 
investigate the relationship between the structural and magnetic phase transition, the dependence 
of the magnetization and the ODW ordering parameters on temperature in the tetragonal and or- 
thorhombic structures is observed. Note that the evolution of sublattice magnetization on tem- 
perature reflects the variation of the spin alignments and the long range order,, and the sublattice 
magnetization vanishes at the critical temperature, characterized by Tsdw'^ while the evolution of 
the 00 parameters associated with the orthorhombic CF splitting reflects the structural phase tran- 
sition, which is characterized by T, . In absence of the orthorhombic distortion, the pure electronic 
interaction leads to Ts < Tsdw, contrast to the experimental results. However, in the presence of 
the orthorhombic distortion, > Tsdw- Therefore, in addition to the e-e interaction, addressing 
the e-1 interaction associated with the orthorhombic distortion clearly is crucial to explain the mag- 
netic and structural phase transitions observed experimentally. Our results support a scenario that 
the magnetic phase transition is driven by the ODW which is mainly induced by the orthorhombic 
lattice distortion in the presence of the intermediate electronic correlation. 

In summary, we have presented the spin and orbital polarizations, as well as SDW and ODW 
ordering, in both the tetragonal and orthorhombic phase, and uncovered the relation between the 
structural and magnetic phase transition. We have shown that the orthorhombic lattice distortion 
is an essential factor to stabilize the SAFM ground state through the formation of ODW ordering. 
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